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Hydrogen- like atoms are of the most important quantum systems solved. Even for 
describing stabilization properties and for testing QED and weak interaction theories 
a great deal can be done at the relativistic atomic physics level (Greiner 1991, Kylstra 
et al 1997, Quiney et al 1997). It is very important then any insight that might 
be given on the properties of hydrogen- like systems. An important tool have been 
the algebraic properties of the set of operators defining the system; these are not only 
connected with the corresponding group and its symmetry algebra but often offer also 
simplified methods for doing some calculations. It is the purpose of this Letter to 
define a new set of operators for the Dirac relativistic hydrogen atom. This comprises 
a non unitary representation of the SU(2) algebra and defines ladder operators for the 
problem. An extra phase is needed to close the algebra but this allows us to solve the 
Dirac hydrogen atom in a neat algebraic way. 

The Dirac Hamiltonian for an hydrogen-like atom is 



H = a-p + Pm , (1) 

r 

where a and (3 are Dirac matrices (Bjorken and Drell 1964), Z is the atomic number, 
r is the relative distance between the electron and the nucleus, m is the mass, and we 
use units such that h = c = 1. 

Taking advantage of the obvious rotational symmetry of H, we can express the 
bound eigenstates of the hydrogen atom as 



'A(M) = if.™''?' I. (2) 



where and 3^ are spinor spherical harmonics of opposite parity and j is the total 
angular momentum (Greiner 1991). It is convenient to define the quantum number 
e such that it equals +1 when / = j -|- 1/2 and it equals —1 when / — j — 1/2 and 
use it instead of parity. The 'big' and 'small' radial components of the bi-spinor (2) 
describing bound states of an hydrogen atom are solutions of the system 



where we have defined the positive definite quantity k := y/m? — E'^, E being the 
energy of the bound state and we have expressed the equations in terms of the 
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dimensionless variable p := kr. For the sake of simplicity we also used C := Ze^, 
Tj := e{j + 1/2), and u := y/{m - E)/{m + E). 

Let us now change to the new variable x defined by p = e^ so that the range of x 
is the open interval (—00, 00) and redefine the radial functions as 

F{p{x)) :=Vm + E [iP_{x) + ^+{x)] , 

, (5) 

G{p{x)) :=Vm - E [•^-(a;) - i^+ix)] . 

With the new functions i(^+{x) and i(^-{x), defining 



/.:=^ + l (6) 
and after some manipulations we get the following system of differential equations 



Mx)=(^rf-C-l^Mx), (7) 
^-{x)=(^Tf-e-\^^-{x), (8) 

for describing an the radial part of an hydrogen atom. As it should become clear 
in what follows, the inclusion of the term 1 /4 in the above equations is necessary to 
close the algebra we purport to construct. Notice that equations (7) and (8) can be 
regarded as an eigensystem with the known eigenvalue 



dx"^ 



+ 2(/x- l)e^ -e' 



2x 



c^:=r?-C-l/4 = i(j + l)-C^ (9) 

as follows from the radial symmetry of the hydrogen atom. 

In order to rewrite the system (7) and (8) making clear the relation with a SU(2) 
algebra let us define the two operators 



^^±^=e^'M|-Te^Ti|: + J), (10) 



dx ' 2 

where we introduced the extra phase ^ besides the 'radial' variable x, and a third 
operator 

^^3 := -i^, (11) 

which depends exclusively on ^; we can alternatively define the two operators Q.\ and 
Q.2 as 
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n, = -{Q+ + Q^), fl2 = - (1^+ - 1^-) . (12) 
z Zi 

The previously defined operators are easily seen to satisfy all the relationships of 
the SU(2) algebra 

[ili, ilj] = ieijk^k, i, j, k=l,2, 3, (13) 

where Einstein summation convention is implied; the algebra can be also expressed in 
terms of the operators fl± 

[^3, n±] = ±n±, and n_] = 2^3, (14) 

which thus play the role of raising and lowering operators. The operator 

:= n-n = nl + nl + nl = -f^- e^- - 2ie^^ - ], (is) 

where O := Oii + 02j + ^^3k, is the Casimir operator since [O^, Oj] = for all i = 1, 2, 3. 

Given the just stated properties, we can choose the operators and O3 and define 
the simultaneous eigenstates VJj{x,^) where u and are, repectively, the eigenvalues 
of fl'^ and O3 — the notation is chosen in analogy with the spherical harmonic 1^™(6', (p) 
case. It must be keep in mind though that we are not here restricted to a compact 
set of parameters given the infinite range of x. Our choice of eigenstates allows us to 
write 

nsV^^{x,o^i^vj;{x,o, n'v^^{x,0-u;vj;{x,c), (16) 

and thus 

Vj;{x,0 = e'^V!:{x), (17) 

where we have tried to use a notation for the a;-function reminiscent of the associated 
Legendre polynomials. 

Now, using equations (7) and (8) we can straightway get il^+{x) — V^{x,^) and 
^|J-{x) = VJ^~^{x, besides, we can easily show that the operators Q± change the 
eigenvalue n to the eigenvalue ± 1, i.e. 



n^v>^{x,0 = c^vjf^\x,a 



(18) 



Letter to the Editor 



5 



where the numbers can be determined from < a;yu|0+0_|a;yU >= C~C^_i and, 
with an appropriate selection of phase, they become = ±-y/ //(^ i 1) ~ ^{^ ~ !)• 
These definitions also establish the conection of our operators with the hydrogen atom 
problem since applying to Vj^{x, ^) or to VJ^~^{x, essentially reproduces equations 
(7) and (8) — as the second of equations (16) clearly exhibits. 

For establishing the Hermiticity (or lack thereof) of the operators introduced, we 
need a scalar product; to this end it suffices to use the following product 



<cj>,i;>= JJ^ P ct>*{^.x)i;{^,x)dx. (19) 

With this interior product, the eigenstates V^{x^^) := jo;,// > form a complete 
orthogonal basis < uj' , ii'\uj^ ji >= Sn^iS^^/. Definition (19) also implies that 0,3 is 
an Hermitian operator, but that Qi = —Q\ and ^2 — that is they are anti- 

Hcrmitian (Martinez-y-Romero et al 1997). Thence is not necessarily positive 
definite; a positive definite operator can be anyway defined as 



fi^-fl = nl-rL\-nl = 2^1 - VL^. (20) 

The action of this operator can be shown to imply that 2ij? > ui, meaning that is 
bounded by below; let us call A this minimum value i.e. A := lA^lmin' With this we 
easily get that ui ~ A(A — 1) and so, since ui = — — 1/4 and we are looking for 
positive A, 



l/^lmin = ^ = «+ 2' (2^) 

where s := +{Tj — C^)^^^- The most important conclusion we can draw from our 
discussion is that A no longer has to be restricted to integer or half-integer values as 
happens necessarily in the standard angular momentum or SU(2) (Hermitian) case. 

This curious and interesting result means that looking for solutions to the Dirac 
hydrogen atom can be regarded also as looking for non-unitary representations of a 
SU(2) algebra labeled by real numbers A — that is, no longer restricted to integer 
or half-integer numbers. In fact, according to equation (21) we have two series of 
eigenvalues depending on whether |U is a positive number or not. In the first case, 
A is the minimum number of an infinite set of positive eigenvalues: // = —A — A;, in 
the second case, —A is the maximum value of the infinite set of negative eigenvalues: 
H = —A — k; in both cases k = 0,1,2,... is a non-negative integer. Notice that 
from the physical point of view the existence of infinite representations of the SU(2) 
algebra makes sense since it is associated with the denumerably infinite set of energy 
eigenvalues of the hydrogen atom. In fact, the energy spectrum of the system is easily 
obtained from (6) and is found to be 
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E 



m 



1 + 



-1/2 



(22) 



in the case of positive eigenvalues we have /U = A + /c = s + fc + 1/2, k = 0, 1, 2, . . ., 
this is precisely the energy spectrum for an Dirac hydrogen atom (Bjorken and Drell 
1964, Greiner 1991). The negative eigenvalues do not led to physically admissible 
eigenstates as we shall see in what follows. 

The ground state of the system follows from the equation f2_|AA >= for the 
positive eigenvalues. The solution of this differential equation can be found to be 



2(A-l/2) 

v/r(2A - 1) 



,(A-l/2)x 



exp(-e^), 



(23) 



where V{y) stands for the Euler-gamma function. As A is the lowest eigenvalue, we 
trivially obtain 



^+{x) =Vx{x), and V-(a;) = 0; (24) 

thus the base state is given by |AA >= e''^^7'^. Using the original variable p, the 
big and the small components for the ground state can be shown to behave as 
F{p) oc \/m + Ep^e'f and G{p) oc ^m — Ep^e'P and, in the negative eigenvalue case, 
the solution behaves as ~ p'^e^ resulting in a divergent behaviour as p — > oo making 
it unsuitable as an eigenfunction of the hydrogen atom; it is obvious that every other 
negative eigenvalue function is also unsuitable. 

The excited states are obtained applying $7+ succesively to |AA >; they end being 
polynomials multiplied by the weight factor p^~^/^e~^ which assures the appropriate 
behaviour of the eigenstates both as p — > as well as p — > oo. A more detailed 
discussion of these new polynomials and graphical representations of them are given 
in a more detailed paper (Martinez- y- Romero et al 1997). 

The main conclusion of this Letter has to do with the infinite-dimensional non- 
unitary representation of SU(2) and where each of the basic operators, excepting Q.^, 
are also non-Hermitian. For example, the matrix elements of are 



< a;/x|ri3|a;//' >= pi5^^>, (25) 

where ^ = ±(A + k), /c = 0, 1, 2, . . ., so its trace vanishes and the determinant of an 
element of the group of the form exp(in3)^ is always 1. The other two operators have 
as their only non-zero matrix elements 

< ujp\ni\u;p ± 1 >= ^^^p(p ± 1) - A(A - 1) (26) 
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and 

< u;fx\n2\u;ii±l >= -^Va*(//± 1) - A(A- 1). (26) 

This means that the trace of both fii and ^2 vanish and that the determinant of group 
elements is one only when the parameter takes imaginary values. We notice also that 
^3 can be given the physical interpretation of producing infinitesimal changes in the 
phase of the state \\u) >= e'^^Pj^(x); that is, it can be associated with the unitary 
operator IJ = e^^-^^ which changes the phase of any state. 

In summary, we have constructed an SU(2) algebra for the relativistic hydrogen in 
the Dirac formulation where the corresponding group is not necessarily compact. We 
must pinpoint that in order to close the algebra we were forced to introduce an extra 
parameter ^ which plays the role of a phase. One of the most noteworthy features 
of the representations reported here is the mixing of a spinorial angular momentum 
character (implying an equally spaced spectrum) with the energy requirements of the 
problem (requiring a differently spaced spectrum); the interplay of these two spectral 
requirements is basically reflected in equations (15), (20), (21) and (22), and in the 
fact that the eigenvalues of the system (7), (8) follows from both the generic radial 
symmetry and the specific features of the interaction. From equation (6) we may also 
notice that in the limit of vanishing interaction, ^ — > 0, the representation collapses 
and, in this case, ji = 1 always. This behaviour is precisely as expected because 
there is no longer any restriction over the eigenvalues and thus the spectrum becomes 
continuous — corresponding to a free Dirac particle. 

The representations of the algebra are labeled by numbers which are neither 
integers nor half-integers as ought to be in the case of the ordinary unitary 
representations. Nevertheless, the operator algebra introduced allows an essentially 
algebraic solutiond of Dirac hydrogen atom which may have various applications (De 
Lange and Raab 1991). It is to be noted also the similitude of our transformed 
equations with the corresponding ones for the case of the Morse potential (Niinez- 
Yepez et al 1995, 1997) and the conections that all our formulation has with systems 
with hidden supersymmetric properties (Benitez et al 1990a,b, Martinez-y-Romero et 
al 1991, Martinez- y- Romero and Salas-Brito 1992, Haag et al 1976) as we will discuss 
at lenght in a forthcoming article. 

This work has benefitted of the comments of R Jauregui, L F Urrutia, A Gonzalez- 
Villanueva and H N Nunez- Yepez and has been partially supported by CONACyT 
(grant 1343P-E9607). ALSB also acknowledges the friendly support of R Micifuz, P 
A Terek, Q Gris, F Sadi, B Kot, B Caro, and U Sasi. 
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